Machian Cosmological Solution in the 
Generalized Scalar- Tensor Theory of 
Gravitation 



A. Miyazaki * 

Department of Physics, Loyola University, New Orleans, LA 70118 

and 

Faculty of Economics, Nagasaki Prefectural University 
Sasebo, Nagasaki 858-8580, Japan 



Abstract 

The Machian cosmological solution satisfying (j) = 0{p/u) is dis- 
cussed for the homogeneous and isotropic universe with a perfect fluid 
(with negative pressure) in the generalized scalar-tensor theory of 
gravitation. We propose cij(</>) — 2) for the coupling function 

in the Machian point of view. The parameter ^ varies in time very 
slowly from ^ = to ^ = 2 because of the physical evolution of matter 
in the universe. When ^ — > 2, the coupling function diverges to — oo 
and the scalar field 4> converges to G'^ . The present mass density is 
precisely predicted if the present time of the universe is given. We 
obtain /jq = 1-6 x 10~^^ 5. cm~^ for to = 1-5 x 10^*^ yr. The universe 
shows the slowly decelerating expansion for the time- varying coupling 
function. 
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In the previous paper we discussed the Machian cosmological solution 
satisfying = 0{p/uj) for the homogeneous and isotropic universe with a 
perfect fluid in the Brans-Dicke theory |0. We found that the gravitational 
constant really approaches to constant when the coefficient 7 of the equation 
of state goes to —1/3 (negative pressure) for the closed model. If we assume 
the present mass density po is identical to the critical density pc, taking 
\uj\ ~ 10^ 1^ into account, we obtain the difference e of the coefficient 7 
from —1/3 (e = 37 + 1) has a value of order 10~^ to support the present 
gravitational constant. These parameters lead to the time- variation of the 



gravitational constant 



G/G 



~ 10 yr ^ which is compatible with the 

recent observational data 

The closed model {k = +1) of this Machian solution is valid for the 
coupling parameter u satisfying u/2 + B < 0, where B = — 3/(^ + 2)(^ — 2) 
and ^ = 1 — 37. The above parameter e ~ 10~^ gives B ~ 10^ which 
means u < —10^ at present. The expansion parameter a{t) is explicitly 
a linear function of t, but its coefficient depends on the parameter ^. As 
the parameter ^ goes to 2 through the quasi-static process, the coefficient 
of the expansion parameter increases very slowly for the constant coupling 
parameter uj. Thus the universe exhibits the slowly accelerating expansion 
at present, which is also compatible with the recent measurements [Q of the 
distances to type la supernovae. 

The above results strongly suggest that the coupling parameter u of the 
Brans-Dicke scalar field should vary in time as the universe expands. A 
tentative conjecture is like uj{t) ~ -8(0- The coupling parameter u is almost 
constant and varies in time very slowly through the quasi-static process of 
the parameter ^. Therefore the time- variation of the coupling parameter 
depends on the physical evolution of matter in the universe. 

In the present paper, we discuss the Machian cosmological solution and 
the behavior of the coupling parameter uj{4>) in the generalized scalar-tensor 
theory of gravitation We consider the homogeneous and isotropic 

universe filled with a perfect fluid with pressure. We require a proper cosmo- 
logical solution satisfles the postulate that the scalar fleld should have the 
asymptotic form cj) = 0{p/ uj) when the coupling parameter uj is large enough 
in order that the scalar fleld converges to zero when |ci;| — > 00. 

The action for the generalized scalar-tensor theory of gravitation 

is described in our sign conventions (c = 1) as 

S = j d^a;v^[-0i?+ 167rL„ - ^g^- cjy^^cp^,] , (1) 
where R is the scalar curvature of the metric g^^,, 4'{x) is the Brans-Dicke 
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scalar field, uj{(j)) is an arbitrary coupling function, and Lm represents the 
Lagrangian for the matter fields. The variation of Eq.(|l|) with respect to g^i, 
and (j) leads to the field equations 



-U(t),^,,-g,un<P), (2) 



87rT + 



(3) 



3 + 2iu{(f)) 

where T^^, is the energy-momentum tensor, T is the contracted energy- 
momentum tensor, and □ denotes the generally- covariant d'Alembertian 
□0 = 0''^^. These field equations satisfy the conservation law of the energy- 
momentum 

T.f = . (4) 
The line element for the Friedmann-Robertson- Walker metric is 

ds^ = -dt" + a^{t)[dx^ + (j\x){de^ + sin^^V)] , (5) 

where 

{sinx for k = +1 {closed space) 
X for k = {flat space) (6) 
sinhx for k = —1 (open space) . 

The energy-momentum tensor for the perfect fluid with pressure p is given 
as 

Tf^U = -PQlJiV - (P + P)Uf,Uy , (7) 

where p is the mass density in comoving coordinates and is the four 
velocity dx^/dr (r is the proper time). The nonvanishing components are 
Too = ~Pi = —pQii (i 7^ 0), and its trace is T = p— 3p for the homogeneous 
and isotropic universe. 

The energy conservation Eq.(^) leads to the equation of continuity 

p + 3-(p + p) = 0. (8) 
a 

Let us suppose the equation of state for the perfect fluid 

p{t) = ^p{t) , (9) 
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where —1^7^ 1/3. Now we are rather interested in the "negative" pres- 
sure. Taking the equation of state into account, we can integrate the equation 
of continuity Eq.(|^) and get 

pit) a"" it) = const, (10) 

where n = 8(7 + 1), which has n = 4, n = 3, and n = 2 for the radiation 
era, the matter-dominated era, and the negative pressure 7 = —1/3 era 
respectively. 

The independent field equations which we need solve simultaneously are 

^ (^2 (A ,</>,87rp Stt (p - 3p) 1 

^la+fcj-^l^l +0 + — - 3 + 2^(0) 0' 



and 



a 

+ 3-C 



Svr (p — 3p) 



du{4>) -2' 



'121 



a 3 + 2cj(0) 

According to the postulate = 0{p/uj), we expect a solution described 

as 

where is an unknown function of t and may generally include u as the 
following form 

m = ^oit) + Oil/iu). (14) 

In the Brans-Dicke theory, we know that $(t) should not include u for the 
perfect fiuid model 1,1. 

The most important thing is to determine an arbitrary coupling function 
uj{(f)). There are many literatures, for example [|l3]"[|Ill- We seek some pos- 
sibilities in the Machian point of view. First, if the coupling function u;(0) 
obeys 

'-^4>' = Sn,Mt) . (15) 



we can adopt the similar procedure for the Machian cosmological solution 
in the Brans-Dicke theory. The coefficient 1] is an arbitrary constant. 

Let us assume that the following relation is satisfied approximately enough 
because the coupling function uj{(f)) is almost constant or large enough: 
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We need check later whether this simphfication is vahd or not for the obtained 
Machian solution. Taking Eq.([T6|) into account, we get from Eqs. ([T^) and 



$ + 3-$ = {^ + ri)p, 



(17) 



where ^ = 1 — 87 or ,^ = 4 — ra. The function $(t) should not include u 
similarly in the Brans-Dicke theory for the Machian solution because the 
right-hand side of Eq.(^3|) does not depend on u. Thus the ratio d/a should 
not also include uj, and so we find for the expansion parameter 

a{t) = A{uj)a{t) , (18) 

where A and a are arbitrary functions of only u and t respectively. 

Taking Eqs. ([T3|) , (|16|) , and (|18]), we obtain from Eq. (|Tl|) after elimination 
of 0by Eq.dllD 



2 



Ap 



3k 



AHuj)a^ 




(3 + r/)p 



(19) 



We are not interested in the flat space {k = 0) here. For the closed and 
the open spaces (fc = ±1), if we require that Eq. (p!9|) is identically satisfied 
for all arbitrary values of u, we find that the coefficient A{uj) must have the 
following form 



^(0) 



+ B 



A^ioo) 

where i? is a constant with no dependence of u. Thus we get 

2 



(20) 



2 



-] +^-k^^ 




(3 + ^)P I u^B 



(21) 

where we introduce a notation j = —1 for uj/2 + B < and j = +1 for 
u/2 + B > 0. Again, by requiring that Eq.(|2T|) is identically satisfied for all 
UJ, we obtain two identities, 



and 





(3 + r7)p 



(22) 



(23) 
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We are convinced of the existence of the following type of solutions on 
the analogy of the previous paper 

$(^) = Cp(^)^^ (24) 

a{t) = bt, (25) 
where coefficients ( and b are constants respectively, and in fact we find 



- (g_2)[g4(?+2)g-2,?] > for kj = -landO^^ <2 
I+M forkj = +land2<^^4:, 



(27) 



(5~2)K2+{r?+2)C-2,y] 



and 



^ r^e^ - (5r/ + 3)e + 

(e - 2) + (r/ + 2)e - 2ry] ^ ^ 

from Eqs.(p!7D, (p2[), and ( p3D successively. 

For the closed space (A; = +1), we require 77 > for ^ < when ^ ^ < 2 
to give the attractive gravitational force {G > 0). On the other hand we 
require 77 < 0, for the coefficient b is real when ^ ^ < 2. So we encounter 
a contradiction. Let us adopt the other alternative that the coefficient rj is 
a linear function of ^ {v = Vi^ + ^o)- We put rjQ = for the finite ( when 
,^ = 0. Thus we replace Eq.(|l5D to 



= 8vrr/iep(t) = Sttt^i (p - 3p) . (29) 

After similar discussions, we obtain as a solution 

C = (r/i + l)/(e-2), (30) 

,2^/ -l/(e-2)(e-ei), forkj = -land0^^<2 

\ l/(e-2)(e-6), forkj = +land2<^^4, ^^'^ 

and 

B = [Viie - 5^ + 3) - 3]/(^ - 2) [(r/i + 1)^ - 2(r/i - 1)] , (32) 

where ^1 = 2(771 - 1)/ (r/i + 1). 

We require —1 < rji < 1 for ( < and the real b when ^ < 2 
(A; = +1). By analyzing the time-dependences of the solution, we get from 
Eq.(ll) 

do; ^ Via3 + 2u) 1 

# (^, + l)(e-2)0- ^ ^ 
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After integration, we obtain the coupling function 



|3 + 2cu(0) I = C(f){t) [2^i?/('7i+i){€-2)] ^ ^34) 

where C is an integral constant. However, we realize that this coupling 
function uj{(f)) is not consistent with Eq.(|T6|) because uj{(j)) varies in time too 
rapidly. 

As the next alternative, we suppose 

d0 ~ 3 + 2u{<P) 

to suppress the rapid time- variation of uj{4>). In this case, the right-hand side 
of Eq.(|T7|) includes the coupling function, and so the function may also 
include as Eq.(|l^. From now on, we restrict our calculations to the first 
order in 1/(3 -|- 2a;). At least for the present time {\uj\ > 10^), this restriction 
would give a good approximation. For the first order in 1/(3 + 2uj), the 
function $(t) does not include uj{(f)). Thus we obtain the Machian solution 
for the first order in 1/(3 + 2uj) again, 

C = - 2) , (36) 

(4 - e)-'^' , forkj = -land0^^<2 
(^2_4)-i/2^ for kj = +land2 <^^A, ^ ' 

and 

i? = -3/(e-2)(e + 2). (38) 

In Eqs. (|33D or (0), the scalar field includes the coupling function c<j(0) 
itself. This expression causes ambiguities, and so let us define another func- 
tion of $ 

uj{(j)) = . (39) 

Taking this expression into account, we estimate Eq.(p5D by the Machian 
solution for the first order in 1/(3 + 2u;) and get 

— = (40) 

which leads to 

w($) = ^ln|$(t)| + lV, (41) 

where W is an integral constant. 

When ^ — >• 2, the time-dependence of p{t) goes to t^^ and so p(t)t^ 
const. However, the coefficient ( diverges to the minus infinity as ^ — 2 in 
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the range of ^ < 2. So the function = C,p{t)t^ also diverges to the 
minus infinity very slowly as ^ 2 in the same range. Thus, if > 0, the 
coupling function vj{^) diverges to the minus infinity when ^ — »• 2 in this 
range. On the other hand, when t — >■ and ^ — *■ 0, <I>(t) and w{^) also 
diverge to the minus infinity. The time-derivative of the scalar field (j){t) is 
described in the present solution as 



0(t) = — ?!L^$(t) 

' 3 + 2ti7(<l>) ^' 



^ 16 7] 



(42) 



3 + 2ti7(<l>) ^ - 2j ' 

and the second term of the middle bracket (~ 1/ln |$(t)|) goes to zero when 
^ — > and ^ — >■ 2. Therefore the assumption Eq.([l6|) is consistent in our 
Machian solution near ^ = and ^ = 2. The coupling function Eq. (^ID 
satisfies the constraint 

^ + B<0. (43) 

which leads to 

r/ln|$(t)| + (e-2)W^>6/(e + 2), (44) 

if 7] > in the range of ^ < 2 (at least near ^ = and ^ = 2, or for the 
appropriate W). The behavior of the scalar field becomes 

when ^ — > 2. The scalar field is approaching to const as the parameter ^ ^ 2 
and its value gradually converges to zero. 
The expansion parameter is expressed as 



a{t) = A{uo)ht 
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ti7($)(e + 2)(e-2)-6. 



1/2 



t (46) 



for the closed space {k = +1) in the range of ^ ^ < 2. Taking ro($) — > — oo 
when t ^ into account, we get a{t) ^ near ^ = 0. By replacing t] to ?7i^, 
we might be able to improve the situation, that is, 

duj{4>) -2 ^ Sttt]! (p - 3p) 
d<P ^ 3 + 2a;(0) ' ^ ' 

and thus 

t^($) = ^ln|$(t)| + iy. (48) 

When we consider the case near ^ = 2, we encounter a severe trouble: the 
coefficient of a{t) in Eq.(^BD converges to zero as $(t) — > —oo when t ^ oo. 
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To satisfy the constraint Eg. (^31) , the couphng function must decrease 

faster than B ~ — 2) when ,^ ^ 2. On the other hand, the couphng 
function must decrease slower than — 2) when ^ — ^> 2 in order to 

give the slowly accelerating expansion. After all, there are no alternatives 
but 

(which would give the almost linear expansion,) for the Machian cosmo logical 
solution. This means 

dzu 

5? = °- <^°) 

The coupling function w{^) does not depend on $. Therefore the solutions 
( P^ ) and ( P5| ) with Eqs.(|5BD, (P7D, and ( PBD are exact for all coupling param- 
eter uj. The constraint Eg. (|^) for all ^ (0 ^ ^ < 2) requires the condition 
T] > 3, which gives uj < —3/2 when ^ = and avoids the singularity. For this 
coefficient r/ > 3, the expansion parameter a{t) gives the extremely slowly 
decelerating expansion except the early stage with = 0. We obtain a(t) = t 
for ^ — > 2 if 7] = 3. 

The scalar field (f) for the coupling function tt7($) = ?7/(,^ — 2)is given as 
the following and we get when ^ 2 

- 2) + 2r7 Tj 

which converges to a definite and finite constant in the limit. If we assume 
?7 = 3, taking to = 1-0 x 10^*^ and the present gravitational constant 
Go = 6.67 X dyn.cim? .g~^ into account, we can estimate the present 
mass density po = 3.5 x 10~^^ (7. cm^^. If we adopt to = 1-5 x 10^° yr, we 
obtain po = 1-6 x 10^^^ g.cm"^, which is very near to the critical density pc ~ 
10~^^ g.cm~'^. As the parameter ^ ^ 2, the coupling function tu($) diverges 
to the minus infinity and the gravitational constant approaches dynamically 
to the constant Goo- 

It should be remarked that the time- variation of the coupling function 
tu($) = ri/{^~2) is derived from that of the parameter ^. However, we have 
not known the details of the physical evolution of matter in the universe yet. 
Probably, our universe started (classically) from the Big Bang with ^ = 
(the radiation era), via the matter-dominated era (^ = 1), and now must be 
approaching to the negative pressure era (the quintessence era, 7 = —1/3, 
^ = 2). According to the recent measurements of the coupling parameter 
(|a;| ~ 10^), we get e = 2 - ^ ~ lO'^ from Eq.(|Up. Over 10^° yr, the state 



of the universe has been varying extremely slowly from ^ = to e ~ 10 . 
This situation would solve the problems existing in the solar system and the 
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evolution of life in the time- varying gravitational constant. At i ~ 5 x 10^ yr, 
the parameter e had already been small enough and the gravitational constant 
did not differ much from the present value. 

There is a discontimiity at ^ = 2 for the Machian cosmological solution, 
and so it is not reasonable that the universe continues to evolve beyond ^ = 2 
to the state ^ = 4. Our universe must be approaching to the state ^ = 2 
for ever {t +oo). We require that the universe is closed {k = +1 ) for 
this range of ,^ (0 ^ ^ < 2) to give the attractive gravitational force. As the 
parameter ^ varies in time extremely slowly (the quasi-static process) in this 
range, we may regard that the parameter ^ is constant when we execute the 
derivative with respect to t. 

It is believed that the generalized scalar-tensor theory of gravitation re- 
duces to general relativity with the same energy-momentum tensor when 
\uj\ -|-oo and u^^duo / d(j) — > 0, but this is not generally true. The present 
Machian cosmological solution has no ordinary correspondences in general 
relativity even when |a;| -\-oo and uj~^duj/d(t) — > 0. The scalar field of 
Machian solutions satisfies the asymptotic form = 0{p/uj) and converges 
generally to zero when |c<j| +oo or p — ^ 0. It should be noted that the 
scalar field 0(t) converges to the finite in the present Machian solution 
when u) — >• —oo and p — >• as t — >• -|-oo, though it satisfies 4> = 0{p/u}). The 
parameter uj is no more arbitrary but is determined by the evolution of the 
universe in the generalized scalar-tensor theory of gravitation. 

When we deal exactly with the time-varying coupling function w{(^) = 
ri/{^ — 2), the expansion parameter a{t) does not show the slowly accelerating 
expansion but rather the extremely slowly decelerating expansion (almost the 
linear expansion). If this is not compatible with the recent observations, we 
might need to introduce furthermore the cosmological constant to the present 
Machian cosmological model. 

We have no criterions to determine the coupling function uj{(j)) in the 
generalized scalar-tensor theory of gravitation. We discussed some examples 
here. They seem to suggest that it is not an arbitrary function of (p but 
another scalar field which is derived from matter itself by another field equa- 
tion. The scalar field (f) was introduced to the Brans-Dicke theory through 
the relation GM/c^R ~ 1. By the similar analogy, the relation u ~ 1/(^ — 2) 
seems to require the existence of another unknown scalar field connected with 
matter. 

The next straightforward problem is to determine the time- variation and 
the range of the parameter ^{t). We need investigate the evolution of a scalar 
field as the dark matter in the universe. 
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